Spectroscopic and magnetic mirages of impurities in nanoscopic systems 
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We present exact results for magnetic impurities in nanoscopic systems with focusing properties. We 
analyze the spectroscopic and magnetic properties of Kondo, intermediate valence and magnetic im- 
purities on a sphere with a metallic surface. Exact calculations show the occurrence of spectroscopic 
and magnetic mirages at the antipodes of the impurity location. Comparison with calculations per- 
formed using effective models validates previous calculations of spectroscopic mirages. Our results 
can be extended to other systems with focusing properties like quantum corrals. 



PACS numbers; 72.15.Qm, 73.22.-f 

Interference of electron wave-functions is a common 
phenomena in condensed matter physics which leads to 
a plethora of effects. Some examples are the formation 
of energy bands in crystals, the Friedel oscillations in the 
screening of charged impurities or the Aharonov-Bohm 
effect in mesoscopic circuits. Up to now, due to techno- 
logical difficulties, it was not possible to design and build 
devises that focus electron wave- functions to artificially 
create images in a metal, in a recent and notorious ex- 
periment, Manoharan et al.EI were able to project the im- 
age of a Kondo impurity to a remote location. The devise 
used in this experiment is an elliptical quantum corral on 
a Cu surface. A Co atom on an open Cu surface behaves 
as a Kondo impuritAjrand generates a characteristic spec- 
troscopic signatureE'H This signature can be observed by 
measuring the tunneling current to a scanning tunneling 
microscope (STM) tip placed at the vicinity of the impu- 
rity. The signal disappears when the tip is moved away 
of the Co impurity a distance of the order of 10 A. When 
a Co atom is placed at one of the foci of the elliptic cor- 
ral, the Kondo signature is not only observed when the 
STM tip is placed close to the impurity but also when it 
is located at the empty focus of the corral indicating the 
coherent refocusing of the electronic structure from one 
focus to the other. 

These experiments triggered a number of theoretical 
works where the eleiutcpnic image projection in quantum 
corrals is analyzed.Q □ All calculations incorporate the 
Kondo effect in a phenomenological way, either by in- 
cluding a Kondo resonance at the Fermi energy or adding 
a phase shift to the conduction electrons. The quantum 
mirages so obtained show that in order to reproduce the 
spectroscopic properties at the empty focus, correlations 
do not play an essential role, they only provide the sig- 
nature of the impurity at the Fermi level where the STM 
experiment probes the local density of electronic states. 
The Kondo nature of the Co impurity is relevant only 
to generate a resonance at the Fermi level. A complete 
many body theory is still lacking and it is of central im- 
portance in order to understand some aspects of the prob- 
lem. There are still many open questions particularly 
concerning the spin dynamics in these systems: While in 
the vicinity of the impurity, conduction electron spins are 



screening the Co spin, are spins at the empty focus fluc- 
tuating coherently with the impurity spin? Moreover, in 
the presence of a magnetic non-Kondo impurity that cou- 
ples ferromagnetically with conduction electrons placed 
at one focus, is there a magnetic mirage at the empty 
focus? 

The aim of this work is to present a many body cal- 
culation of impurities in nanoscopic systems with focus- 
ing properties and analyze local density of states at the 
Fermi level and spin correlations. The calculation is per- 
formed by exact diagonalization in small clusters using 
the Lanczos algorithmaE3. The elliptic corral is not the 
only geometry with focusing properties, other examples 
are two mutually facing parabolas at a Cu surface or a 
thin metallic film covering a rugby ball. In the latter 
case, an impurity placed at one tip generates a mirage at 
the other. A particular case of this geometry is a sphere 
with a metallic surface where any point can be defined as 
a focus and the mirage is formed at the opposite point: 
given a point that defines a pole, all meridians cross at 
the other pole being the "optical " length of all of them 
the same. The main difference between these geometries 
is the degeneracy of the one electron states, however a 
deep understanding of the refocusing properties in one of 
them will shed some light into all the others. 

We analyze the case of the sphere that is the simplest 
case and is appropriate for exact diagonalization. Con- 
duction electrons are confined at the surface of a sphere 
of radius _R, the electron coordinates are the angles (0, 1^9) 
and a Kondo impudty is placed at 9 ~ Q. The Anderson 
Hamiltonian readald: 

l,m.a CF 

+ II^(cLc?, + 4qo.) (1) 

where c|^^ creates an electron with quantum numbers 
m,tT corresponding to the total angular momentum, the 
z-component of the angular momentum and the spin re- 
spectively and with energy e; = l{l + 1) /2mR^ , the op- 
erator d^. creates an electron with spin a at the impurity 
orbital with energy Ed and Coulomb repulsion U. The 
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last term describes the hybridization with Vi = VYio{0, 0) 
where Yim{0,ip) is the spherical harmonic. For an im- 
purity at the pole {9 ~ 0), only the m — states are 
hybridized with the impurity orbital. In the summation 
over I we take a cutoff /max to avoid including wave func- 
tions oscillating with a characteristic wave length shorter 
than a typical interatomic distance a, ((a/i?)Zmax — !)■ 
In what follows we refer to the impurity site (9 = 0) as 
the occupied focus and to the point 6* = tt as the empty 
focus. Energies are taken in units of /2mE? and we 
consider a very large on-site Coulomb repulsion U . 
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FIG. 1. Spectral density for Ed ~ —20 and two values oi 
the hybridization: V/4tt = 1 ((a), (c), (e)) and V/47r = E 
((b), (d), (f)). From top to bottom: impurity spectral density 
conduction electrons spectral density for different angles in ar 
expanded scale around sf and conduction electrons spectra 
energy difference with and without impurity. In (e) and (f J 
curves have been shifted to coincide at ef- The Fermi energj 
is indicated with a vertical dotted line. 

We first present exact numerical results for the loca 
density of states for a system with /,„ax — 10 and 14 elec- 
trons. In Fig.(l) the impurity spectral density is showr 
for different values of the parameters (a Lorentzian widtl: 
of 0.5 is taken). A small hybridization V corresponds tc 
the Kondo regime with an average number of electrons 
in the impurity orbital < rid > close to one. For the pa- 
rameters of Fig. (la) we obtain < >= 0.961 and the 
impurity spectral density Ai(ijj) shows a large peak at 
an energy close to the impurity energy Ej^ and a number 
of small peaks at the band energies. Within the band, 
the dominant peak coincides with the Fermi energy sp 



and represents the Kondo resonance. A large hybridiza- 
tion {V > \Eii\) reduces the number of electrons a the 
impurity orbital and for the parameters of Fig (lb) we 
obtain < Ud >= 0.705 corresponding to an intermediate 
valence (IV) regime. In this case Ai{Lu) shows a large 
amplitude at the band energies and does not present a 
dominant Kondo peak. The conduction electron spec- 
tral density A{ui, 9) is shown also in Fig. (1) for different 
values of the angle 9. At the impurity focus, the con- 
duction electron spectral density presents a pseudogap, 
both for the Kondo and IV regimes (Figs (Ic) and (Id)). 
As 9 increases the pseudogap is filled and at the equa- 
tor [9 = 7r/2) the local density of states resembles the 
unperturbed density of states. As 9 approaches tt, the 
empty focus, a clear mirage of the impurity is obtained 
both for small and large V . In this simple geometry, the 
states that fill the pseudogap for 9 ^ Q,tt correspond 
to the states with m ^ 0. For a better comparison, in 
Figs, (le) and (If) we present the difference 6A{u!,9) — 
A{u!,9) — A^{u!,9) of the conduction electrons spectral 
density with and without the impurity. This is a typical 
result when the Fermi energy coincides with one of the 
non-interacting one-electron levels. When the Fermi en- 
ergy lies between two conduction energy levels i.e. when 
there is a closed shell in the sphere, the Kondo peak is 
not observed and the conduction electron spectral den- 
sity shows a small dependence with 9. These results qual- 
itatively reproduce the experimental observations in the 
quantum corral where the mirage is observed only for 
system sizes such that a confined one-electron state lies 
close to the Fremi energyJll 
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FIG. 2. Conduction electron spectral densities and spectral 
density difference for an effective resonant level at ef. A 
small system as in Fig.(l) ((a), (c)) and a large system with 
/max = 50 and 220 electrons ((b), (d)). 

The low energy spectrum (uj « ep) obtaind for the 
Kondo impurity can be mimicked by representing the 
Kondo resonance as an uncorrelated resonant level at the 
Fermi energy, which is precisely the spirit of some pre- 
vious calculations where the iiripurity propagator is not 
fully self-consistently evaluatecH. A simple calculation 
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using Hamiltonian (1) with = ep , U = and a small 
effective hybridization V gives the results of Fig. (2): 
for a small system the spectral density for oj k, ep qual- 
itatively reproduce the exact results of Figure (1). Due 
to the simple structure of the one electron wave func- 
tions on a sphere it is clear that, although the structure 
of the difference 5A{u), 9) on the two foci for energies 
close to the Fermi energy is different, at the Fermi energy 
6A{eF0) = SA{eF, tt) independently on the radius of the 
sphere provided that the Fermi level coincides with one 
of the one-electron levels. An alternative approximate 
procedure is to express the conduction-electron propaga- 
tors in terms of the impurity propagator Gi{io) and then 
use for Gi{uj) a Lorentzian peak at the Fermi energy with 
a width given by the Kondo temperature Tk-o However, 
in this procedure Gi{uj) does not account for the struc- 
ture of the discrete energy levels of the nanoscopic system 
and the exact results obtained for an isolated small sys- 
tem are not reproduced. This approach is justified only 
if there is a large broadening of the one-electron levels 
of the nanoscopic system, a limit in which we expect a 
smaller and more diffuse mirage. 

These results validate some previous results where the 
Kondo resonance is not fully self consistently calculated 
and the many body effects are taken into account in a 
rather phenomenological way. In addition our results 
show that also for IV impurities, with a large hybridiza- 
tion or a small Ed, the mirage in the spectral density 
is still obtained. Another very interesting feature con- 
cerns the spin correlations in these systems. We start by 
presenting the static spin correlation < S^mp • cr(^, > 
where Simp is the impurity spin operator, a{0, ip) is the 
conduction electron spin operator at coordinates (0, ip) 
and the brackets indicate the expectation value at the 
ground state. For a small hybridization the results pre- 
sented in Figure (3a) clearly show an antiferromagnetic 
correlation for small resulting from the screening of 
the impurity spin by the conduction electrons. As 9 
grows, the spin-spin correlation decreases, however, at 
the empty focus it increases again giving rise to a mag- 
netic mirage. At the empty focus, the conduction elec- 
tron spin is fluctuating coherently with the impurity spin 
as if the Kondo impurity were in its neighborhood. When 
the Fermi energy coincides with one of the one electron 
levels, the Kondo effect takes place and the conduction 
electron spin at the empty focus is always antiferromag- 
netically correlated with the impurity spin. For the case 
of a closed shell, the spin-spin correlations are weaker 
and at the empty focus the conduction electron spin is 
ferromagnetically correlated with the impurity spin. For 
a large hybridization charge fluctuations are important 
and although the local {9 — 0) spin-spin correlations are 
large the magnetic mirage is not well defined (dashed line 
in Fig. (3a)). 
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FIG. 3. Spin-spin correlations as a function of 9. (a) An- 
derson Hamiltonian (1); (b) Kondo Hamiltonian (2) with fer- 
romagnetic and antiferromagnetic coupling J; (c) Ising-type 
interaction normalized to the value at 6 — for a small 
(dashed) and large (continuous line) system. Squares are the 
exact results obtained in the full Hamiltonian (2). 

In the Kondo regime the Anderson Hamiltonian can be 
approximated by the Kondo Hamiltonian mliich describes 
correctly the spin dynamics of the systemll3: 
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a(0,0) 



(2) 



Lm,a 



here J > is the antiferromagnetic coupling between 
the impurity and the conduction electron spins. For 
J < 0, this Hamiltonian describes the effect of a mag- 
netic nooJCondo impurity with a ferromagnetic (FM) sd- 
couplingtll. We have also calculated the spin correlations 
for a magnetic non-Kondo impurity with FM coupling to 
the conduction electrons. The results obtained for J > 
and J < are compared in Fig. (3b). The spin correla- 
tion at short distances (small 9) are larger for the Kondo 
impurity. This is due to the singular nature of the Kondo 
scattering which is not obtained for a ferromagnetic sd- 
coupling. The magnetic mirage is proportionally more 
pronounced for the latter case where perturbations or 
simple approximations can be used to estimate the spin- 
spin correlations. For the J < we have calculated the 
correlations with an Ising-type interaction J^Sf^^a^ (0, 0) 
and, as shown in Fig. (3c), the results quantitatively re- 
produce the exact results obtained with the full Hamil- 
tonian This allows us to use, for the case J < , the 
Ising-type interaction to estimate the spin correlations 
in much larger systems and in more complicated geome- 
tries. The spin correlations obtained for a large sphere 
are also shown in Fig. (3c). The magnetic mirage is quite 
pronounced and despite the fact that the two foci can be 
even further than 50 interatomic distances, the spin cor- 
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relation for electrons at the empty focus are as large as 
the correlations at distances of the order of one or two 
interatomic distances from the occupied focus. 

A systematic analysis of the behavior of the magnetic 
mirage can be made by changing the sphere radius R 
and keeping the electron density constant. As the ra- 
dius increases, the characteristic energy difference SE be- 
tween two consecutive one-electron levels decreases and 
for some particular values of R the Fermi level coincides 
with one of the one-electron states. We analyzed these 
cases in which there is an open (non-filled) shell in the 
sphere. The results obtained for different radii are pre- 
sented in Fig. (4) , the inset shows the amplitude of the 
spin-spin correlation between the impurity spin and the 
conduction electron spin at the empty focus for fixed den- 
sity (same as in Fig. 1) and Jz = 0.25 as a function of 
R. At short distances {9 « 0) the spin-spin correlation 
is sensitive to both, the system size and the cutoff /,nax 
and as increases it becomes independent of /max- At 
the empty focus two different regimes can be observed: 
For small systems, the level spacing SE is smaller than 
Jz and the impurity polarizes the electron at the Fermi 
level; the magnetization profile is essentially given by the 
one-electron wave function \Yip^m=oi9)\'^ where Ip is the 
Fermi angular momentum. Since these wave functions 
at 9 = TT increase as Ip increases, the magnetic mirage 
increases with R for small systems. For large systems, 
SE « Jz, the impurity produces a mixing of states with 
different I that interfere destructively and the magnetic 
mirage tends to cero as R goes to infinity. 
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FIG. 4. Spin-spin correlations vs 6 for the Ising limit for 
different R with fixed Jz and electron density (open shell). 
Inset: Spin correlation at S = tt as a function of R for the 
same Jz 

As a consequence of the presence of the magnetic mi- 
rage, we can infer that if an additional impurity is placed 
at the empty focus (i.e. when there are two impurities in 
the system, one at each focus), they will strongly interact. 
In this case the impurity-impurity coupling is ferromag- 
netic if the Fermi level lies close to one of the energy 
levels of the nanoscopic system and is antiferromagnetic 



and much weaker if the Fermi energy lies between two 
one-electron levels. Results obtained for the magnetic 
mirage in the elliptic quantum corral will be presented 
elsewhere. 

In summary, we have presented here an exact many 
body calculation for impurities in nanoscopic systems 
with focusing properties. The spectroscopic character- 
istics of the Kondo impurity, which are always observed 
in its neighborhood for open shell systems, are repro- 
duced at the opposite point on the sphere giving rise to 
a spectroscopic mirage of the impurity. Our results also 
show that impurities with a large hybridization, which 
induces large charge fluctuations, also give rise to spec- 
troscopic mirages. This may be important since ab-initio 
calculations of Co on noble metal surfaces indicate-that 
charge fluctuations are relevant for these systemal^.We 
compared the exact results for the Kondo impurity with 
those obtained with an effective model which includes a 
resonant state at the Fermi level. Although there are 
some quantitative differences between the two models, 
the exact results are qualitatively reproduced by the sim- 
ple one-body approach. We have also calculated the ex- 
act spin-spin correlations in small systems. The results 
show the existence of magnetic mirages: at the empty fo- 
cus conduction electron spins are fluctuating coherently 
with the impurity spin. These magnetic mirages are also 
obtained for non-Kondo magnetic impurities which cou- 
ple ferromagnetically with conduction electron spins. For 
this case the magnetic mirage is more pronounced than 
for the Kondo impurity. 

The results obtained here for a sphere can be extended 
to other systems with focusing properties such as quan- 
tum corrals. In particular, in addition to a mirage in the 
localized density of states, a magnetic mirage should be 
observed in such systems in the presence of a Kondo (or 
sd) impurity. Results for the elliptical corral will be pre- 
sented elsewhere. 
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